Critical frequency for vortex nucleation in Bose-Fermi mixtures in optical lattices 



in 
o 
o 

> 

o 



5h 



O 

o 



> 
o 



in 
o 

-I— » 



i 

-a 
c 

o 
o 



X 
S3 



M. Guilleumas, M. Centelles, M. Barranco, R. Mayol, and M. Pi 
Departament d'Estructura i Constituents de la Materia, Facultat de Fisica, 
Universitat de Barcelona, Diagonal 647, 08028 Barcelona, Spain 

We investigate within mean-field theory the influence of a one- dimensional optical lattice and of 
trapped degenerate fermions on the critical rotational frequency for vortex line creation in a Bose- 
Einstein condensate. We consider laser intensities of the lattice such that quantum coherence across 
the condensate is ensured. We find a sizeable decrease of the thermodynamic critical frequency 
for vortex nucleation with increasing applied laser strength and suggest suitable parameters for 
experimental observation. Since 87 Rb- 40 K mixtures may undergo collapse, we analyze the related 
question of how the optical lattice affects the mechanical stability of the system. 
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An essential feature of a rotating superfluid is the even- 
tual appearance of quantized vortices |lj . For this reason, 
nucleation of vortices in dilute gases in rotating traps has 
been the subject of intensive experimental and theoreti- 
cal research (see e.g. Refs. [ilHIillHIilHIljIIlllil and 
references quoted therein). 

Motivated by recent advances in the manipulation of 
rotating Bose-Einstein condensates (BEC) Oil 0] and 
of condensates loaded into optical lattices |13(. and in 
anticipation of future experiments, we study a slowly 
rotating BEC in a one-dimensional (ID) optical lattice 
[l4|. Rotating ultracold Bose gases in optical lattices 
are suitable systems to study quantum phenomena. In 
particular, the excitations, stability and dynamics of a 
vortex line in the superfluid regime have been recently 
addressed 0. Concurrently, the physics of degenerate 
Bose-Fermi mixtures in optical lattices is attracting con- 
spicuous attention. Theoretical progress is being made in 
the basic understanding of the system's quantum phase 
diagram and the superfluid to Mott-insulator transition 
|16|. and experimentalists have already prepared a de- 
generate mixture of Rb and K atoms in a tight optical 
lattice [l7j . 

The aim of the present work is twofold. Firstly, to 
address the formation of vortices in a coherent array of 
Bose-Einstein condensates in a ID optical lattice, and to 
determine the dependence of the thermodynamic critical 
angular velocity £l c for vortex formation [2j,|3j] on the laser 
intensity. Secondly, to investigate the effect of trapped 
fermions on £l c . Attractive Bose-Fermi mixtures may 
experience collapse and, thus, we also discuss the impact 
of the optical lattice on the mechanical stability of the 
mixture. 

We consider a zero temperature mixture made of a 
87 Rb Bose condensate (B) and a degenerate 40 K Fermi 
gas (F). They are confined by the axially-symmetric ex- 
ternal potentials of a harmonic magnetic trap and of a 
stationary periodic optical lattice modulated along the 
z-axis. The lattice is produced by a far detuned laser, 
which hinders the possibility of spontaneuous scattering 
and yields practically equal Rb and K optical potentials. 
The resulting potential for each kind of atom q = B,F 
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where m q is the atomic mass. The radial and axial fre- 
quencies of the harmonic trap are taken from a recent 
experiment 13]: ujb± = 2itvb±_ with Vbx = 74 Hz and 
ijjBz = %'kvbz with vbz = 11 Hz for 87 Rb, while those for 
40 K are a factor [ttib /mp) 1 / 2 — 1-47 larger. The optical 
potential is determined by its period d = 2.7 '/im [13J and 
depth Vq — sEr, where Er — h 2 n 2 /2md 2 = h x 80 Hz is 
the recoil energy and s is a dimensionless parameter. 

For 3 x 10 5 atoms of 87 Rb confined in the harmonic 
trap, the condensate is cigar-shaped with a Thomas- 
Fermi length Ltf = 84 /jm, and a radius -Rtf = 6 /jm. 
The superimposed optical potential splits the BEC over 
Ltf/c? ~ 30 wells. We shall be concerned with a range 
of laser intensities that ensures full coherence of the con- 
densate across the whole system |l5|. The criterion for 
the Mott transition for the considered parameters 
leads to an estimate of a maximum lattice depth for the 
superfluid regime of Vb ~ 150 Er ~ h x 12 kHz |lj| , when 
quantum tunneling between consecutive wells is still suf- 
ficient to retain coherence. 

Within our mean-field approach, neglecting p-wave in- 
teractions, the energy density functional that describes 
the boson-fermion mixture at zero temperature with a 
singly quantized vortex in the condensate along the z- 
axis has the form fl9ll2CI 
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where ns — \^\ 2 is the condensate density, is 
the fermion density, and Tp is the fcrmion kinetic en- 
ergy density written as a function of nj? and its gradi- 
ents [Thomas- Fermi- Weizsacker (TFW) approximation]. 
The boson-boson and boson-fermion coupling constants 
gsB and gsF are written in terms of the s-wave scat- 
tering lengths as and cibf as gsB — 4-7T agli 2 /ms 
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FIG. 1: (Color online) Contour plots of the 87 Rb condensate 
density in the xz-plane. (a) Condensate in the pure harmonic 
trap, (b) Same as (a) with a quantized vortex along the z- 
axis. (c) Same as (a) under an optical lattice with Vo/h = 
2.5 kHz. (d) Same as (c) with a quantized vortex along the 
z-axis. Distances are in units of the oscillator length obj_ = 
iJh/(mBUB±) = 1-25/ttm. The density scale from to 1 at 
the top of the figure is normalized to the highest density in 
each panel [these are 229 in (a), 212 in (b), 394 in (c), and 
376 a~ 3 ± in (d)]. 



and gBF = 47r a,BFh 2 /ttibf , respectively, with ttibf = 
2mBniF I (tub + mp). 

The variation of £ with respect to tib and np , un- 
der the constraint of given number of bosons Nb and 
fermions Np, yields two coupled Euler-Lagrange equa- 
tions: a Gross-Pitaevskii equation for bosons with a term 
describing the boson-fermion interaction, and a TFW 
equation for fermions 0, |2(| , which is increasingly valid 
the larger Np is. Starting from randomly sampled den- 
sities, we solve these equations as indicated in Ref. 12 1| . 
We use the adapted version of the virial theorem |2(| to 
check the numerical convergence. 

We first consider a pure BEC, with Na = 3 x 10 5 atoms 
of 87 Rb as in the experiment of Ref. For the boson- 
boson scattering leng th we use the value as = 98.98 ao 
(lao = 0.529 A) j 7 ^. In Fig. ^ we show contour plots 
of the condensate density in the xz-plane, ns(x,0, z) = 
\^(x, 0, z)\ 2 , for different situations. In (a) the conden- 
sate is confined only by the harmonic trap. In (b) the 
cigar-shaped condensate of (a) hosts a quantized vortex 
line along the z-axis. In (c) a ID optical lattice with 
Vo/h = 2.5 kHz splits the cigar-shaped condensate (a) 
into an array of multiple disk- like coherent condensates. 
As shown in Ref. [2^ , certain aspects of the macroscopic 
properties and low-energy dynamics of a magnetically 
trapped BEC in a tight optical lattice can be understood 




FIG. 2: (Color online) Density profile of the condensate along 
the z-axis ns(0,0,z) = |^ (0, 0, z)\ 2 , in units of ffljj*^, for dif- 
ferent laser intensities. The coordinate z is in units of obi. 
Inset: column density 723(2;), in units of 10 3 x ag 2 ± , for the 
same Vq values. The column density of the condensate when 
it hosts a vortex line is also plotted. The coordinate x is in 
units of a B ±- 



in terms of a renormalized interaction coupling constant 
g* BB > §bb and of an effective mass m* B > tub along the 
direction of the periodic optical potential. The increase 
of the radial size of the sample that one observes in panel 
(c) with respect to panel (a) is due to the increased re- 
pulsive effect of the boson-boson interaction when the 
system feels the optical lattice. It originates from the 
fact that the optical confinement produces a local com- 
pression of the gas inside each well 0, |2^| . The axial 
size of the condensate also increases appreciably due to 
the combined effect of the repulsive interactions and the 
redistribution of atoms inside the potential wells. 

In panel (d) of Fig. ^ the optical lattice is superim- 
posed to the condensate with the vortex state (b). In 
a bulk superfluid, the size of the vortex core is of the 
order of the healing length £ = [87mBa] -1 / 2 , where ub 
is the bulk density. This expression holds for inhomoge- 
neous superfluids taking for ns the local boson density 
in the absence of vortices y|. Thus, the vortex core size 
is larger for the outer sites of the split condensate, where 
the density is smaller. 

Figure [21 shows the evolution of the density pro- 
file of the condensate along the z-axis, nsiQ, 0, z) = 
1^(0, 0, z)| 2 , with the laser intensity. In the presence of 
a shallow ID optical lattice, the condensate profile starts 
developing small oscillations around the profile at Vq = 
which follow the periodic optical potential, but its spatial 
extension is nearly not affected. When the laser intensity 
increases the system nearly splits into separate conden- 
sates. The effect of the local compression inside each po- 
tential well of the lattice, as well as the increase of the ax- 
ial size of the system, can be clearly seen in Fig. EI In the 
inset we plot the column density ub(x) = J dzns{x, 0, z) 
for the condensate in the harmonic trap alone and when 
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FIG. 3: (Color online) Critical angular velocity Q c for nucle- 
ation of a singly quantized vortex line as a function of the 
laser intensity Vq. 



the optical lattice is added. We also display n b (x) for the 
condensate hosting a vortex line, in this case the quan- 
tized circulation around the vortex line pushes the atoms 
away from it. 

The thermodynamic critical angular velocity £l c for nu- 
cleating a singly quantized vortex is obtained by sub- 
tracting from the vortex state energy E\ in the rotating 
frame the ground-state energy E , i.e., f2 c = (Ei/Nb — 
Eq/Nb)/^, 0, 0, and it provides a lower bound to the 
critical angular velocity In Fig. 03 we plot fl c as a 
function of the laser intensity for Nb = 3 x 10 5 and for 
Nb = 5 x 10 4 . For a fixed number of condensate atoms, 
O c decreases when Vq increases, in agreement with the 
fact that the radial size of the system becomes larger 
in the presence of the optical lattice and the associated 
reduction of atoms along the symmetry axis (see the de- 
crease of the column density ns at the origin in Fig. [21 
with increasing Vb). In turn, the increase of the inertia 
of the condensate along the direction of the laser beam 
due to the larger effective mass m* B |23j also contributes 
to the reduction of the critical angular velocity. 

Assuming that the relative effect remains of the same 
order in experiment, Fig.[3]predicts a sizeable 24% reduc- 
tion of fl c already with a laser strength Vo/h = 5 kHz for 
both of the Nb values. The effect is enhanced the shal- 
lower the magnetic trap is. A less elongated magnetic 
trap would favor the appearence of a vortex at a slower 
rotation. For example, the values of O c (Vb = 0) shown in 
Fig. would be decreased by a 0.6 factor in a spherical 
trap with Vbz = vb±_ = 74 Hz. The relative reduction of 
£7 C caused by Vo/h = 5 kHz with respect to Vb = would 
still be of 20-23%. 

Figure|31suggests an experiment to test the dependence 
of the critical frequency for vortex nucleation on the in- 
tensity of the optical lattice, since it shows that, for a ro- 
tating condensate with a given Nb, a vortex line should 
nucleate at angular frequencies lower than f2 c (Vb = 0) if 
a co-rotating ID optical lattice is superimposed, whereas 
otherwise it would not. As the thermodynamic f2 c un- 
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FIG. 4: (Color online) Contour plots of the condensate (B) 
and fermionic (F) densities in the xz-plane for the 87 Rb- 
40 K mixture, in the combined harmonic and optical lattice 
trap. The condensate hosts a quantized vortex. In (a) and 
(b) V /h = 0.5 kHz, whereas in (c) and (d) Vo/h = 2.5 kHz. 
Distances are in units of as±- The density scale from to 1 
at the top of the figure is normalized to the highest density 
in each panel [these are 312 in (a), 25 in (b), 457 in (c), and 
37 a- 3 ± in (d)]. 



derestimates the actual critical frequency for vortex cre- 
ation, dynamical calculations are needed for a precise 
quantitative estimate of this effect. 

With the purpose of assessing the effects of a trapped 
fermion cloud on the rotating condensate in an opti- 
cal lattice, we consider next a 87 Rb- 40 K mixture. It is 
characterized by a large attractive Bose-Fermi interac- 
tion that assists the sympathetic cooling of the fermionic 
species down to the degenerate regime. There are vari- 
ous values of the interspecies scattering length as f pub- 
lished, depending on the typ e of experiment and tech- 
niques used El HI 123 . A recent measurement 
of asF at JILA [27j has established a value asF = 
—250 ± 30 ao, which we adopt here because of its reduced 
experimental uncertainties. Vortex states in 87 Rb- 40 K 
mixtures confine d by a harmonic potential have been re- 
cently addressed [2C| (although asF — — 395 ao was 
employed). When a ID optical lattice is switched on, the 
atomic species experience the trapping potential (|T|). 

We have carried out calculations for a mixture with 
N B = 3 x 10 5 and N F = 1.5 x 10 5 atoms, and the 
same parameters as in Fig. ^ We assume that the Fermi 
component is in the normal — nonsupcrfiuid — but quan- 
tum degenerate phase, and consider that it is in a sta- 
tionary state. This situation could be achieved experi- 
mentally by waiting long enough after the generation of 
the vortex in the condensate to let the drag force be- 
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tween bosons and fermions to dissipate. We present in 
Fig. |1| contour plots of the boson n B (x, 0, z) and fermion 
np(x, 0, z) densities when the condensate hosts a vortex 
line, for two laser strengths: Vo/h = 0.5 (a,b) and 2.5 
kHz (c,d). The condensate and the fermionic cloud are 
modulated by the regular pattern of the optical lattice. 
The large mutual attraction between fermions and bosons 
makes the effective interaction between bosons less re- 
pulsive than in the pure BEC. Furthermore, the boson 
atoms induce an effective attraction between fermions 
counteracting the Fermi pressure. As a result, the den- 
sity of both species increases in the overlapping region, 
as if they were more strongly confined by the external 
potentials, and the condensate becomes more compact 
in space (compare Figs. 0t and^l). Still, the effect of 
the Pauli exclusion principle is notorious in Fig. 0] where 
the 40 K cloud is seen to extend to larger distances than 
the 87 Rb atoms, both axially and radially. 

The presence of a vortex along the z-axis in the con- 
densate component is apparent in panels (a) and (c) of 
Fig-EJas the boson density vanishes on the vortex line, ft 
is interesting to note that the large attractive cibf leads 
to a visible depletion of the fermionic density on the z- 
axis as well, reminiscent of the bosonic vortex core. The 
effect should be directly observable experimentally, ft is 
more evident in (d) where bosons and fermions are more 
tightly confined by the optical lattice. 

In the present mixture we find that the presence of 
the quantum degenerate fermionic cloud does not change 
much the value of Q c compared to the pure BEC, nor 
its pattern against the optical lattice strength. Indeed, 
for all the Vo values considered here the critical fre- 
quency for vortex appearence is raised by ~ 10%, due 
to the enhancement of the condensate density at the 
core caused by the fermion atoms. The impact of the 
fermion cloud on O c is magnified in a deeper magnetic 
trap. For instance, if the same mixture was set in a trap 
with vbz = 95 Hz and vb_l = 640 Hz, that preserves 
the aspect ratio vbz/vbi^ = 0.15 used in Fig. 0] the crit- 
ical frequency f2 c (Vb = 0) would increase by 40% with 
respect to the Np — case. In this squeezing magnetic 
trap, however, the effect of Vq on fi c is negligible (< 2%). 

A peculiar feature of the 87 Rb- 40 K system, stemming 
from the interplay of the moderately repulsive as and 
the strongly attractive clbf, is the existence of a mechan- 
ical stability limit beyond which the system cannot sus- 
tain more fermions, for at high densities the strong inter- 
species attraction may overcome the Pauli pressure and 
drive the ultracold gas mixture into collapse. The phe- 
nomenon has been observed experimentally |25j. Thus, 
we next address the stability of the 87 Rb- 40 K system sub- 
ject to an optical lattice. In fact, collapse is currently the 



focus of experimental attention because of the implica- 
tions the instability has for constraining clbf |25l l26l |27| . 
and also because there are prospects that a radially- 
squeezed 87 Rb- 40 K mixture at a density close to collapse 
must be able to form stable bright soliton trains |2j| . 

Keeping Nb fixed at 3 x 10 5 atoms as in the previ- 
ous discussions, with clbf = — 250ao we find that the 
boson-fermion mixture confined in the magnetic trap of 
Fig. 0] would be stable up to virtually arbitrary fermion 
numbers, consistently with the findings of Ref. [Tflj . In- 
creasing u>b±_ prompts the occurrence of collapse. In a 
trap with vbz = 95 Hz and vbi_ = 640 Hz, the 3 x 10 5 
rubidium atoms are able to retain up to a maximum of 
7V™ ax = 2.6 x 10 5 fermions, when the system collapses. 
(The effect on Q, c of Np = 1.5 X 10 5 and V for this trap 
has been discussed above.) Now, a decrease of 7V™ ax can 
be seen when a ID optical lattice is added, caused by the 
effective local compression of the atoms inside the opti- 
cal trap. Indeed, superimposing a ID optical lattice with 
d = 2.7/im and Vo/h = 5 kHz to this harmonic trap, the 
87 Rb- 40 K mixture with Nb — 3 x 10 5 can sustain up to 
A™ ax = 7.5 x 10 4 fermions only, i.e., less than 30% of 
the Vo = value. The effect appears to be accessible 
to verification under present experimental conditions. It 
opens the possibility to study the collapse of a mixture for 
trapped fermion and boson numbers considerably smaller 
than in a pure harmonic trap. 

In conclusion, we have determined the thermodynamic 
critical frequency for nucleating a quantized vortex in an 
array of coherent condensates within the Gross-Pitaevskii 
theory. We have studied a Bose- Fermi mixture trapped in 
a ID optical lattice to examine the effect on the value of 
fl c of a fermionic cloud in the quantum degenerate phase. 
87 Rb- 40 K mixtures may collapse and, in this regard, it 
has been shown that a ID optical lattice can be an effi- 
cient degree of freedom to tune the onset of instability. 
Our analysis of the relative variation of Q c calculated by 
thermodynamic arguments is intended to provide a useful 
information for future experiments. A quantitative de- 
termination of f2 c requires dynamical calculations. Work 
in this line is being undertaken. 
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